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Abstract 

TlCuCla has a singlet ground state with the excitation gap A = 0.65 meV. 

The magnetic excitations in TlCuCls have been investigated by means of 

neutron inelastic scattering experiments. The constant- Q energy scan profiles 

were collected in the a* — c* plane. A well-defined single magnetic excitation 

mode was observed. The dispersion relations along four different directions 

were determined. The lowest excitation occurs at Q = {h, 0, 1) with integer 

h and odd /, as observed in KCuCls. A cluster series expansion up to the 

sixth order was applied to analyze the dispersion relations, and the individual 

exchange interactions were evaluated. It was demonstrated that TlCuCls is 

a strongly coupled spin-dimer system. 

PACS numbers: 75.10.Jm, 75.30.Et, 75.40. Gb, 78.70.Nx 
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I. INTRODUCTION 

The magnetic excitations in quantum spin systems, which have a singlet ground state 
with an excitation gap, are a new subject in magnetism. Recently, neutron inelastic scat- 
tering experiments have been actively performed to investigate the magnetic excitations in 
various coupled spin-dimer systems, e.g., (VO)2P207 [|ll, Cu(N03)2 ■ ID2C) 0, BaCuSi206 
0, SrCu2(B03)2 [3 and KCuCls A variety of dispersion relations for the propaga- 

tion of the excited triplet have been observed in these systems, and the true nature of the 
exchange networks, which was not expected from the crystal structures, has been revealed 
in (VO)2P207 and KCuClg. 

Since large single crystals of KCuCls can be obtained, its magnetic excitations have been 
extensively investigated by neutron inelastic scattering [p|-p^. The dispersion relations were 



first analyzed using the effective dimer approximation, in which the exchange interactions 



between individual dimers are reduced to an effective interaction between dimers [p,p,10,ll 



Later, the cluster series expansion was applied by Miiller and Mikeska to describe the 
dispersion relation in order to evaluate the individual exchange interactions. From these 
analyses, the exchange network in KCuCls was elucidated. Consequently, KCuCls has been 
characterized as a weakly and three-dimensionally coupled spin-dimer system. 

In this study, we investigate the magnetic excitations in TlCuCls, which is isostructural 
with KCuCla |13|. This compound has a monoclinic structure (space group P2i/c) [jl4[. 
The crystal structure is composed of planar dimers of CU2CI6, which are stacked on top 
of one another to form infinite double chains parallel to the crystallographic a-axis. These 
double chains are located at the corners and center of the unit cell in the b — c plane, and are 
separated by T1+ ions. Figure 1 shows the projection of Cu^"*" ions with spin-| on the a — c 
plane and the exchange interactions between Cu^"^ ions. For the notation of the exchange 
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interaction, see reference IT^ . 



The magnetic ground state of TlCuCls is a spin singlet with an excitation gap, as observed 



in KCuCls ||T3|. The magnitude of the spin gap A in TlCuCls was evaluated from the critical 



field of the magnetization curve |T5|,|T6[ and the excitation energy of the direct ESR transition 
as A = 0.65 meV. The lattice parameters, the critical fields and the saturation fields for 
KCuCls and TlCuCls are listed in Table I. The crystal lattice of TlCuCls is compressed along 
the a-axis and enlarged in the 6 — c plane as compared with KCuCls. Thus, substituting 
Tl"*" for K"*" produces uniaxial stress along the a-axis. 

Although the crystal structures of TlCuCls and KCuCls are the same, there is a sig- 
nificant difference between their magnetic properties, i.e., the spin gap for TlCuCls is one- 
quarter of that for KCuCls, while the saturation field for TlCuCls is about twice as large as 



that for KCuCls |T5|JT8|1 . This suggests that the interdimer interactions in TlCuCls are much 
stronger than those in KCuCls. However, the details of the individual exchange interactions 
in TlCuCls have not been clarified so far. 

Recently the field-induced magnetic phase transition was observed by magnetization 
and specific heat measurements |TB|,|TP[. It was demonstrated that the nature of the phase 



transition can be described in terms of the Bose-Einstein condensation (BEC) of the triplet 
excitations (magnons) |2^. In the magnon BEC theory, the magnons around the lowest 
excitation are relevant to the phase transition. For these reasons, we have investigated the 
magnetic excitations in TlCuCls. 



Recently Cavadini et al. |2T]] have also independently investigated the magnetic excita- 
tions in TlCuCls by means of neutron inelastic scattering. They measured mainly the dis- 
persions parallel to the principal axes a*, h* and c*, and evaluated the intradimer exchange 
interaction and the effective interactions between dimers, which are given by the linear 
combinations of the individual exchange interactions as shown in section III. Although the 
present measurements are confined in the a*c*-plane, we measured the dispersion relations 
not only parallel to the a* and c*-axes, but also for two diagonal directions (/i, 0, 2h + 1) 
and (/i, 0, —2h + 1.4), which are roughly perpendicular to each other. As shown in the next 



section, the magnetic excitation is most dispersive along {h, 0, 2h+l), while it is less disper- 
sive along [h, 0, —2h + 1.4). The {h, 0, 2h + 1) direction is parallel to the cleavage (1, 0, —2) 
plane, in which the hole orbitals of Cu^"^ spread. The dispersion relations along {h, 0, 2h+ 1) 
and {h, 0, ~2h + 1.4) are essential, because the interdimer exchange interactions cannot be 
uniquely determined without them, i.e., there is another set of exchange parameters which 
can fit the dispersion relations principal axes a*, b* and c*. We analyzed the dispersion 
relations using a cluster series expansion to sixth order to evaluate the individual interdimer 
exchange interactions, which cannot be determined within the framework of the effective 



dimer approximation [^. The effective interactions calculated with these individual inter- 
dimer exchange interactions are significantly different from those obtained by effective dimer 
approximation. 

The arrangement of this paper is as follows. In the next section, experimental details 
are described, and the experimental results are presented. In section III, the experimental 
results are supplemented by a theoretical analysis based on dispersion curves calculated from 
a cluster series expansion, and the individual exchange interactions are determined. Section 
IV is devoted to conclusions. 

II. EXPERIMENTS AND RESULTS 

TlCuCla single crystals were grown from a melt by the Bridgman method. The details of 
sample preparation have been reported in reference [|16|. The TlCuCla crystal cleaves easily 



parallel to the (0, 1, 0) plane, in which lie the a*- and c*-axes. 

Neutron inelastic scattering was performed using the ISSP-PONTA spectrometer in- 
stalled at JRR-3M, Tokai. The constant-/cf mode was taken with a fixed final neutron 
energy Ef of 14.8 meV. In order to gain intensity, collimations were set as open - monochro- 
mator - 80' - sample - 80' - analyzer - 80' - detector. The energy resolution was about 2 meV 
because of loose collimations. A pyrolytic graphite filter was placed after the sample to sup- 
press the higher order contaminations. We used a sample with a volume of approximately 
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2.5 cm^. The sample was mounted in an ILL-type orange cryostat with its a*- and c*-axes 
in the scattering plane. The crystallographic parameters were determined as a* = 1.6059 
1/A, c* = 0.71513 1/A and cos /3*=0. 0967 at helium temperatures. 

In the previous neutron inelastic scattering measurements for TlCuCls [Q, three exci- 
tations were observed in the energy range E < 15 meV for the scans along {h, 0, 0), {h, 0, 1), 
(1,0,/), (1.5,0,/) and {h,0,2h - 1) with l</i<1.5 and 0</<l. From the temperature de- 
pendence of the excitation spectra, we concluded that the lowest dispersive excitation is of 
magnetic origin. Since the phonon excitation has been observed at almost the same energy 
in the isostructural KCuCla ||23|, we inferred that the highest dispersionless excitation can 
be attributed to the phonon excitation. However, the origin of the second excitation was 
unclear. 

We first investigated the temperature variation of these excitations for Q = (1.5, 0, 0). At 
T = 1.5 K, three excitations were observed at ~ 3, 7 and 12 meV. At room temperature, 
which is much higher than the temperatures corresponding to the excitation energies, the 
peak intensity of the lowest excitation decreases to the background level, while the intensities 
of the two higher excitations appear to be almost independent of temperature. This indicates 
that the two higher excitations are not intrinsic magnetic excitations. If the cause of the 
excitations is phonons, their intensities should increase with increasing temperature. Hence, 
it is too early to conclude that the two higher excitations are attributable only to phonons. 

In the present study, we focused on the lowest magnetic excitation, and investigated the 
precise dispersion relations along (/i, 0,1), (0,0,/), {h,0,2h + 1) and {h,0, —2h + 1.4) by 
constant-^ energy scanning for < h < 0.5 and 1 < / < 2 (see Fig. 2). The present scan 
area is closer to the origin of the reciprocal space than the previous one, so that the magnetic 
excitations can gain intensity due to the magnetic form factor. The dispersion relations along 
diagonal directions {h, 0, 2h+l) and {h, 0, — 2/1+1.4) are necessary to determine the exchange 
network uniquely within the data for the a* — c* scattering plane. 

Figure 3 shows the scan profiles for Q = {h, 0, 1), (0, 0, /), {h, 0, 2h+l) and {h, 0, -2h+lA) 
measured at T = 1.5 K. A well-defined single excitation can be observed in almost all scans. 
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The scan profiles were fitted with a Gaussian function to evaluate the excitation energy, as 
shown by the solid lines in Fig. 3. The horizontal bars in Fig. 3 denote the calculated 
resolution widths. Almost all peaks have widths equal to the resolution limit. 

Figure 4 shows the constant-Q scans for Q = (0.1,0, 1.2) measured at T = 1.5 K and 
80 K. A single excitation is clearly observed aX E = 6.3 meV and T = 1.5 K. At T = 80 
K, the excitation spectrum broadens out to the background level. This indicates that the 
excitation can be attributed to a magnetic origin. 

The dispersion relation uj{Q) obtained for the a* — c* scattering plane is summarized in 
Fig. 5. It is evident that the lowest excitation occurs a.t Q = (0,0, 1). However, since the 
excitation energy is lower than 1 meV, we could not determine the excitation energy due 
to the incoherent scattering and the low energy resolution. In Fig. 5, we substituted the 



gap energy A = 0.65 meV, which was evaluated from the previous magnetization |15,16 



and ESR measurements |T^, for the excitation energy at Q = (0,0, 1). Recently Cavadini 
et al. |21|] investigated the magnetic excitations in TlCuCla by neutron inelastic scattering 
experiments. Their experimental results along the a*-and c*-directions are in agreement 
with our results. 

Based on the present results, together with the previous ones, it is evident that the 
periodicity of the magnetic excitation in TlCuCls is the same as that of the nuclear reciprocal 
lattice along the a*-axis, but doubled along the c*-axis, as observed in KCuCla . Hence, 
it is deduced that the lowest excitation occurs aX Q = {h, 0, /) with integer h and odd / in 
the a* — c* plane. This is consistent with the results of recent neutron elastic scattering 
experiments in magnetic fields [|2^. When a magnetic field H is applied in the present 



system, the single excitation splits into three excitations, since the excitation should be a 
triplet excitation. The lowest excitation energy decreases with increasing magnetic field, 
and finally becomes zero at the critical field = A/g^-Q. For H > the system can 
undergo magnetic ordering due to three-dimensional (3D) interactions. Such field-induced 
magnetic ordering has actually been observed by neutron elastic scattering for H \\ b and 
H > He ^ 5.5 T Magnetic Bragg reflections were observed at Q = {h, 0, /) with integer 
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h and odd / in the a* — c* plane, and are equivalent to those for the lowest excitation at zero 
field. 

There is a sharp contrast between the dispersion curves for Q = {h, 0, 2h + 1) and 
[h, 0, ~2h + 1.4) which are roughly perpendicular to each other. The magnetic excitation is 
most dispersive along {h, 0, 2h + 1), which is parallel to the (1, 0, —2) cleavage plane, and is 
less dispersive for Q perpendicular to it. The dispersion curve for Q = {h, 0, 2h + 1) has a 
local minimum dX h = 0.25. This dispersion behavior is similar to that in KCuCla, which 
implies that the principal exchange pathways in TlCuCls should be the same as those in 
KCuCls. However, the dispersion range in TlCuCls (0.65 meV < u; < 7.3 meV) is much 
larger than that in KCuCla (2.7 meV < oj < 5.0 meV). This suggests that the exchange 
interactions in TlCuCls are much stronger than those in KCuCla. The details of the exchange 
interactions are evaluated in the next section. 



III. ANALYSIS AND DISCUSSION 

The main point in the analysis of the experimental results is that the magnetic interac- 
tions in TlCuCls are dominated by the exchange between the two spins forming the planar 
CU2CI6 dimer (intradimer exchange), whereas exchange interactions between other spin pairs 
(interdimer exchange) can be considered as weaker. The strong antiferromagnetic intradimer 
interaction is the origin of the singlet ground state. 

Since notable anisotropy effects have not been observed in the static measurements of 
TlCuCls |ll6| , p^ , the full magnetic interactions will be described by the spin-| Heisenberg 



model 



n = Y. (5. ■ 5,) . (1) 



For the exchange interactions between spins, we will use the notation given in ref. []I2 
the main intradimer exchange is denoted as J . The exchange interaction per bond between 
spins in dimers separated by a lattice vector la + mb + nc is denoted as the exchange energy 
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Jimn for pairs of spins at equivalent positions in their respective dimer and as J/^„ for spins 
at inequivalent positions. Finite values of interdimer interactions J... will be considered for 
the following exchange paths only: (Imn) = (100), (200), (Oi|) and (201). In the 

following we present an analysis of the lowest elementary triplet excitation based on this 
picture, in three steps. 

(i) The simplest picture is to assume unperturbed propagation of the excitation of one 
dimer from the singlet to triplet state, neglecting the excitation of further triplets during 
propagation. This leads to the basic dispersion 

u^iQ) ^ J + sJ2\Q) , (2) 
sJ2\Q) = [J(foo)cos(27r/i) + J(foo)cos(47r/i) + J(foi)Cos{27r(2/i + /)}] 

±2[J(fi i)Cos{7r(2/i + l)}cos{'Kk) + J^^iifos{'Kk)cos{'Kl) . (3) 

At this level, only certain combinations of the exchange interactions enter the dispersion 
law: these are known as effective dimer interactions: 




The factor in front of the individual exchange constants J and J' results from the number 
of identical exchange paths for the same lattice vector between dimers; it is one of 1 (identical 
exchange on two legs or on two diagonals) , | (exchange on one diagonal only whereas there 
is no exchange path for the other diagonal) or (the exchange path does not contribute). 
Dispersion relation (2) is valid in the limit Jimn, J'lmn ^ J ■ 

Since there are two different dimers per chemical unit cell, the dispersion law has two 
branches, distinguished by ± in eq.(2). However, the structure factor at the zeroth order is 
of the form 
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S{Q, uj±) ~ ^sin ^ ± sin 



QR2V 



(5) 



where Ri = 0.47a+0. 106 + 0. 31c and R2 = 0.47a- 0.106+ 0.31c denote the spin separations 
in the CU2CI6 dimers located at the corner and at the center of the unit cell in the b — c 
plane. Under the present experimental condition, i.e., Q in the a* — c* plane, the uj^{Q) 
branch gives the only nonvanishing contribution. Hence, we can assume that the observed 
single excitation corresponds to the uj^{Q) branch. 

(ii) The simple approach (i) gives a qualitatively correct picture, which confirms the 
assumption of one dominating intradimer exchange interaction, quantitatively. However, it 
is not sufficient. The standard way to improve on this has been to treat the intermediate 



excitation of two or more triplets in an RPA-like approximation ||6|,|25[|. This approximation 
continues to treat dimers as units and results in the dispersion law 



uJ±iQ) = ^J^ + 2J5Jl\Q). (6) 

This result depends only on the effective interactions Jf^^^ ^^e dispersion law does not dis- 
tinguish between the contributions to the interdimer interactions from the different exchange 
paths. 

The measured dispersions were fitted to eq. (6) using all the effective exchange constants 
given in eq. (4). The results are presented in Table II. 

(iii) A systematic approach to improve on the one triplet dispersion law of eq. (2), is 
to expand the energy of the elementary triplet excitation order by order in the interdimer 



exchange interactions. The method has been described in ||T2[, and it has been found that, 
from the second order, new terms which are not present in the expansion of the square 
root of eq. (6) appear. To determine individual exchange couplings, we have implemented 
the cluster expansion technique |T^J2^J^ to allow us to calculate one magnon dispersion 
relations up to the sixth order. For tractability, only exchange paths along (100), (201) and 



were taken into account, whereas the small (as evident from approach (ii)) interactions 



along (200) and (0^3) were disregarded. Due to the complex interaction structure, it was 



not possible to take into account only topologically different clusters, and we had to consider 
all embeddings, i.e., all the different ways of setting a given cluster on the lattice, up to 
the desired order. This leads to an exponentially growing number of cluster from order to 
order. In the present case, we perform calculations for the sixth order treating 18084 colored 
clusters and 673826 embeddings. 

The series expansion gives the dispersion relation in the form 

u{Q)/J=l+ E C^^;,:;;!^,fl«f cos[7rK/i + n,0]cos(7rn,A;), (7) 

pi>0,...,P5>0 {nh,nk,ni) j=l 

where a^'s are the individual exchange interactions in units of J: 

O^l = •^(201)A "2 = J{100)/J, "3 = -^(lOO)A «4 = J(i^^)/J, «5 = i)A (8) 

Leading terms in can be combined to an expansion of a square root, as discussed in 
|l^ . We have calculated the coefficients C^l'nfn, explicitly up to the fourth order, i.e., 



J2%iPj < 4; these coefficients are available up on request. Up to the sixth order, we 
are able to obtain the dispersion numerically for a given parameter set In order to 

demonstrate the convergence of the cluster expansion we show in Fig. 6 the results for 
increasing orders in the dimer expansion for the direction Q = {h, 0, 1). 

For the analysis of the data, we proceeded in two steps. We first performed a least-squares 
fit to the measured data using the analytically known dispersion to the fourth order with the 
intradimer interaction J = 5.68 meV, which was obtained in (ii). Guided by the result of 
this fit, we performed calculations to the sixth order. The parameters for the best sixth order 
calculation are given in Table III. The solid lines in Fig. 5 indicate the calculated results. 
The experimental dispersion curves can be reproduced well by the present calculation. 

The results for the effective exchange constants are well defined and imply, in particular, 
a significant increase in the value of i„ when compared to the RPA-like approach (ii) 
above. The individual exchange constants are well defined for (Imn) = (100) and (201). 
For (Imn) = (l|^) the individual exchange constants are determined less reliably. A fit of 
comparable accuracy is obtained for ^(iii) = 0.23 meV, ^[ii i) = —1-37 meV and minor 
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changes in the remaining parameters. Similar to the situation in KCuCls, the neighboring 
dimers couple magnetically along the chain and in the (1, 0, —2) plane. The most important 
interdimer interaction is the diagonal ^(201) interaction, which is about half of the intradimer 
interaction. 



IV. CONCLUSIONS 

We have presented the results of neutron inelastic scattering for the spin gap system 
TlCuCls. Well-defined magnetic excitation spectra were observed in the a* — c* plane. The 
dispersion relations of the magnetic excitations in TlCuCls were determined as shown in Fig. 
5, and were analyzed by the cluster series expansion. The individual exchange interactions, 
which could not be obtained from the effective dimer approximation, were evaluated, as 
shown in Table III. TlCuCls is a strongly coupled spin-dimer system, in contrast to KCuCls 
which is characterized as a weakly coupled spin-dimer system. The analysis of the individual 
exchange constants shows that the ladder system in TlCuCls, similar to KCuCls, is much 
closer to an alternating spin chain than has been previously believed. 
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FIGURES 




FIG. 1. Projection of Cu ions with spin- 2 on the a — c plane and the exchange interactions. 
The double chains located at the corner and the center of the chemical unit cell in the b — c plane 
are represented by solid and dashed lines, respectively. The shaded area is the chemical unit cell 
in the a — c plane. 



14 




-0.5 



FIG. 2. Scanning directions for Q along (a) (/i,0, 1), (b) (0,0,/), (c) (/i,0,2/i + 1) and (d) 
(/i,0, -2/i + 1.4). 
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FIG. 3. Profiles of the constant-Q energy scans in TlCuCls for Q along (a) {h, 0, 1), (b) (0, 0, 1), 
(c) {h, 0, 2h + 1) and (d) {h, 0, ~2h + 1.4) with < /i < 0.5 and 1 < i < 2. The solid lines are fit 
using a Gaussian function. The horizontal error bars indicate the calculated resolution widths. 
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FIG. 4. Constant-Q energy scans in TlCuCla at Q = (0.1, 0, 1.2) for T = 1.5 K and 80 K. The 
solid line for T = 1.5 K is a Gaussian fit. 
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FIG. 5. Dispersion relations uj{Q) in TlCuCla for Q along (a) (/i,0,l), (b) (0,0,/), (c) 
{h, 0, 2h + 1) and (d) {h, 0, —2h + 1.4). Solid lines are the dispersion curves calculated by cluster 
series expansion to the sixth order using the exchange constants in Table III. 
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FIG. 6. The convergence of the cluster expansion with increasing orders in the dimer expansion 
for Q = {h, 0, 1). Dotted, dashed and sohd Unes denote the results for the second, fourth and sixth 
orders with the exchange constants in Table III. Closed circles are experimental data. 
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TABLES 

TABLE L The lattice constants a, b, c and (3 at room temperature, the critical fields 
He = A/gfiB and the saturation fields Hs for KCuCla and TICUCI3. 



KCuCls 



ref. 



TlCuCla 



ref. 



a (A) 
b{A) 

c(A) 

p 

{g/2)H, (T) 

{g^Hs (T) 



4.029 
13.785 
8.736 
97.33° 
23.1 
54.5 












3.982 
14.144 
8.890 
96.32° 

5.7 
~ 100 



m 
m 



m 
m 

m 



TABLE II. The intradimer interaction J and the effective interdimer interactions Jff^ ^ in 

(Imn) 

KCuCls and TlCuCls. All energies are in units of meV. 

J^ff [meV] KCuClg TlCuCls 

J 4.34 5.68 

Jffoo) -0.21 -0.46 

J(foo) 0.03 0.05 

j;fii. 0.28 0.49 
'■■'-22) 

-0.003 -0.06 

y^2 2' 

Jfl,,) -0.45 -1.53 
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TABLE III. Interdimer exchange interactions in TlCuCla determined by the cluster series ex- 
pansion with the dimer interaction J = 5.68 meV. All energies are in units of meV. AF and F 

denote antiferromagnetic and ferromagnetic exchange interactions, respectively. 

J(ioo) = 0.34 (AF), J(',oo) = 1-70 (AF) J^f^^) = -0.51 

J(iii)= 0.91 (AF), j; -0.57 (F) ^ J,fi i, = 0.74 

•^(201) = 2.56 (AF) J(fo,) = -1.28 
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